Heavy hybrid mesons in the QCD sum rule 
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We study the spectra of the hybrid mesons containing one heavy quark (qQg) within the framework 
of QCD sum rules in the heavy quark Umit. The derived sum rules are stable with the variation 
of the Borel parameter within their corresponding working ranges. The extracted binding energy 
for the heavy hybrid doublets H{S) and M{T) is almost degenerate. We also calculate the pionic 
couplings between these heavy hybrid and the conventional heavy meson doublets using the light- 
cone QCD sum rule method. The extracted coupling constants are rather small as a whole. With 
these couplings we make a rough estimate of the partial widths of these pionic decay channels. 

PACS numbers: 12.39.Mk, 12.38.Lg, 12.39.Hg 

Keywords: Heavy hybrid meson, QCD sum rule, Heavy quark effective theory 



I. INTRODUCTION 



Hadron states which can not be accommodated in the conventional quark model have attracted much interest over 
the past few decades, partly due to the great success of quark model in the classification of hadrons and the calculation 
of hadronic parameters. Theoretically, the existence of these unconventional hadrons may be allowed by Quantum 
II ' Chromodynamics (QCD), the widely accepted fundamental theory of the strong interaction. 

O |, These unconventional hadrons include multi-quark states {qqqq, qqqqq, ■■■)■, glueballs (gg, ggg, ■ ■ ■), and hybrids 
' (995)- Some of them are totally "exotic" , namely their J^*^ quantum numbers are excluded by the conventional quark 
I . model. A straightforward analysis of J^'~^ reveals that 0~~, 0^ , 1 ^, 2"' , • • • are the so called exotic ones. Several 
mesons with exotic J^'~^ = 1 ^, e.g. 7ri(1400) l], 7ri(1600) |2i], have been reported in recent years. These states are 
I ] usually interpreted to be hybrid mesons. The 1 ^ states have been studied in the framework of QCD sum rules in 
^ ■ several works, including their masses [sj and decay properties Q. They were also investigated extensively in other 
\ theoretical schemes, such as Lattice QCD, the flux tube model, and AdS/QCD etc. 

If the above mentioned hybrid mesons exist, there should also be hybrid mesons containing one heavy quark {qQg) 
and heavy hybrid quarkonium {QQg), although the J^*^ quantum number of the former is not exotic. The hybrid 
mesons containing one heavy quark and heavy hybrid quarkonium have been studied in Q . The masses of the hybrid 
' quarkoniums were calculated in the heavy quark limit [oj. The masses and the pionic couplings to conventional 
heavy mesons of the hybrids containing one heavy quark were studied in Ref. [7|]. In the present work, we study 
the hybrid mesons containing one heavy quark in the framework of heavy quark effective theory (HQET) 8], in 
which the expansion is performed in terms of l/niQ, where Q is the heavy quark involved. At the leading order of 
1/toq, the HQET Lagrangian respects the heavy quark flavor-spin symmetry, therefore heavy hadrons form a series 
of degenerate doublets. The two members in a doublet share the same quantum number ji, the angular momentum 
\^ , of the light components. The two ji = ^ S'-wave conventional heavy mesons form a doublet (0~, 1~) denoted as H 
^ ■ and the ji = |/f P-wave doublets (0+, l+)/(l+, 2+) are denoted as S/T. We denote the ji = |/f -D-wave doubtlets 
(1^, 2^)/(2^ , 3^) as M/N. As far as the heavy hybrid containing one heavy quark are concerned, the two ji ~ ^ 
doublets with P = + and P = — are denoted as S'^ and H'^, respectively. Similarly, the two ji — | doublets with 
P = + and P = — are denoted as and M'^, respectively. 

We calculate the binding energy and decay constants of these hybrid doublets using Shifman-Vainshtcin-Zakharov 
(SVZ) sum rules After performing the operator product expansion (OPE) of the T product of two interpolating 
currents, we obtain sum rules which relate the binding energy and decay constants of corresponding hybrid mesons 
to expressions containing vacuum condensates parameterizing the QCD nonperturbative effect. The nonperturbative 
method used to calculate the pionic couplings between heavy hybrid mesons and conventional heavy mesons is the 
light-cone QCD sum rules (LCQSR) 10]. Now the OPE of the T product of two interpolating currents sandwiched 
between the vacuum and an hadronic state is performed near the light-cone rather than at a small distance as in SVZ 
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sum rules, and the QCD nonperturbative effects are included in the light-cone distribution amplitudes of the pion 
state. 

The paper is organized as follows. We construct the interpolating currents for the doublets D^{D = H / S/T /M) in 
Sec. ini We derive the sum rules for the binding energy and decay constants for these doublets in Sec. IIIIl The sum 
rules for their pionic couplings to the doublets H and S are derived in Sec. IIVI The last section is a short summary. 
The light cone distribution amplitudes of the pion which are employed in the present calculation are collected in the 
appendix. 



II. INTERPOLATING CURRENTS 



We adopt the following interpolating currents for the doublets i?'^ and Af 



-hyigs-ffat ■ Gq , 



3Gf 7^ + t-f^at ■ G 



J 



m!} 



Gaip OL-2 



q, 



GOL20 Ctl 



9, 



(1) 



where Gap — GJ^^A"/2 and h^{x) — e"^^'"'^'^^^Q{x) is the heavy quark field with 4- velocity v. The subscript t is 
used to denote the transverseness of the corresponding Lorentz tensors to v, namely 



7t = 7 



9t 
0-t 



«/3 



Gf = G"^ - G°'''vy - G^^^v^v" . 
The overlapping amplitudes between the above currents and the corresponding hybrids are defined as 

{o\jH.mH^{v)) = f„., 

(0|J«.(0)|i7f(T;,A))=/^.7?^.KA), 
(0| J^j. iO)\Mi\v, A)> = Uj.t^Ij, {v, A) , 
{0\J'^TmMi{v, A)) = fM-vTi!^Hv, A) , 



(2) 



(3) 



where ri[v, A) is the polarization tensor of the corresponding heavy hybrid. Apparently these polarization tensors are 
traceless, symmetric to their Lorentz index and transversal to v: 77aQ2 - Qre''^" = 0- Furthermore, summations on A give 
the following projection operators: 



^7?"(^;,A)r7^(^;,A) = -5f , 
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The angular momentum of the light components ji of the heavy hybrids h that couple to the current j\°!h must 



be 3/2. To illuminate this point, we first notice that the 4-velocity of the light quark q of h should be w, namely 

fq = q. This leads to ^i?" = -i?" and 7^7?" = 0, where R°' = G^^-fp + §7^ crj ■ G q. In other words, is a 

Rarita-Schwinger spinor describing a spin 3/2 particle when the 4-velocity of q is v. 

Now it is straightforward to give the interpolating currents for the doublets S'^ and T'' by adding 75 to the currents 
in Eq. 
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J^t = Kg sib 
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The corresponding overlapping amplitudes and projection operators can be defined similarly to Eq. (|18p and (U), 
respectively. 



III. BINDING ENERGY 



To derive the sum rules for the binding energy for the doublets H'^ and M'*, we consider the following correlation 
functions: 

d^xe^'^-^{o\T{jH.{x)4Ao)}\o) = n^.(^) , 

d'xe^'^-{0\T{J^^ix)4i{0)m = -gfu^.iu;) , 
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where uj = 2k ■ v. The correlation functions for doublets 5'' and are similar to that of i/'' and M^, respectively. 

At the quark level(cj ^ 0), it is convenient to calculate the above correlators in coordinate space by OPE at a short 
distance x — > 0. The fourier transformation to the momentum space is straightforward after employing the heavy 
quark propagator in the infinity quark mass limit mg — > oo: 



{0\T{K{x)h40)}\0) = 
The quark propagator used in the OPE of n(w) is 

i A" 



1 + ^ 



dt6{x — vt) . 
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{Q\T{q{x)mm = 
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We consider the contributions of the condensates with dimension not greater than seven in our calculation. The 
contribution involving the quark gluon mixed condensates appears as asiqgsC ■ Gq). It is typically suppressed by a 
factor mQ/{l6T^) compared to the quark condensate, where ttiq = (^gsf • Gq)/{qq) = 0.8 GeV^ and T ^ 1 GeV. 
We thus ignore this contribution. The terms involving the 6-dimensional four quark condensate appear as a'j.{qq)'^ . 
It is of high order in as so it can also be omitted safely. The Feynman diagrams corresponding to the quark-level 
calculation are presented in Fig. [TJ The Fock-Schwinger gauge x^A^^{x) = adopted in our calculation, together with 
the heavy quark propagator in Eq. (jT)), leads to the leading order Lagrangian of HQET £o — hyiv ■ Dhy — hyiv ■ dhy. 
This indicates the decoupling of the heavy quark from gauge field in the heavy quark limit and greatly simplifies our 
calculation. 

The OPE results for n(a;) read 



n^h (w) : 

and similarly 
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(c) (d) (e) 

FIG. 1: The Feynman diagrams for n(w). The double solid line denotes the propagator of heavy quark Q. 



where (GG) = (g^G^^G^/s), (GGG) = {glr^^G^i^G^i^^G''^^), and we used the following formulas in our calculation: 



(gsG'^pG^s) = [gajgi3s - gasg^] {GG) , 



{gsf "(^Ij.uGaffG'p^) = —(GGG) [g^iaga^gisp + gixpgapgav + gaagnpg^f} + gpvgixag0,j 



—giipgacrgpv — g^agapgup — gavg^pgea — g^pgnagvu] 



(10) 



It turned out that n^h(w) = n^(i(w), nj(^h(a;) — \l}^h{uj), Wgh{oj) = Ugjjh{uj), Urph{uj) = Urph{uj), and we 
denote them by 11^-^(0;), njvffe(a;), Ush{uj), and Ilxh{uj), respectively. This implies that the binding energy and 
the overlapping amplitudes of the two members of a doublet are degenerate, which is dictated by the heavy quark 
flavor-spin symmetry. 

For n(w) we have the following dispersion relation: 



n(u;) 



oo 
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With the phenomenological spectral density Pph = f S{s — 2A) H , we can rewrite the above equation to be 
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(11) 
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where A = rrih — mq is the binding energy of the heavy hybrid h containing a heavy quark Q, Pope(s) is the spectral 
density obtained by OPE at the quark level. By performing the Borel transformation 



Blifiu)] = lim 



i-co) 



n+l 



(13) 



aj = -nT 



which is used to suppress the continuum contribution, we obtain the following sum rule with the continuum contri- 
bution subtracted: 



Jo 



(14) 



It is convenient to obtain the spectral density Pope by performing a second Borel transformation B^^^j,[f{T)] to the 
right hand side of the above equation: 



pHh{s) 



38407r3 



327r2 
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FIG. 2: The sum rule for K^h with continuum 
threshold sq = 5-1, 5.3, 5.5 GeV and the working 
interval 0.8 < T < 1.2 GeV. 



FIG. 3: The sum rule for h.gh with continuum 
threshold sq = 5.8, 6.0, 6.2 GeV and the working 
interval 0.8 < T < 1.2 GeV. 
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The binding energy and overlapping amplitudes can now be expressed as 



A 



/;° sp{s)e-'/^ds 
2/;V(s)e-«/^ds 
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In our numerical analysis, we use = 47r/(ll - 2/3n/)ln(so/2/AQCD)2 with = 4 and Aqcd = 220 MeV. 
The mass of up and down quark is ignored. The quark condensate and gluon condensate adopt the standard values 
{qq) = -(0.225 GeV)^, (GG) = 0.038 GeV*. There are several values for the triple gluon condensate (GGG) = 
0.045 GeV"^^], 0.06 - 0.1 GeV^fll], 0.4 GeV^0. The uncertainty caused by this difference is within 6% for A^h 
and Kgh. So we will fix {GGG) to be 0.045 GeV^ in the following analysis. Notice from Eq. ^ and ^ that the only 
difference between the OPE for and or and is the triple gluon condensate. We can then conclude that 
the binding energy of and are almost degenerate, so is the case of and . 

From the requirement that the contribution of terms in the OPE is at least three times larger than that of the next 
term, except for (GGG), we get the lower limit of T denoted by T,„i„. If we require that T^ax > Tmin + 0.4 GeV and 
the pole contribution is at least 20% of the whole sum rule, the lower limit of continuum threshold sq is determined. 

The sum rules for A^h and hgh are plotted in Fig. [2] and Fig. [3l respectively. Finally, we list the extracted the 
binding energy and the overlapping amplitudes in Table ID 
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TABLE I: The values of A, /, and their corresponding sq. 



IV. PIONIC COUPLINGS 



To derive the sum rules for the pionic couplings of these heavy hybrids to conventional heavy mesons, we need the 
following interpolating currents of conventional heavy meson doublets H and S: 
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4r = (17) 

The overlapping amplitudes between the above currents and the corresponding heavy mesons are 

{Q\JHM\Hoiv)) ^ fHo, 
(0|J^,(0)|i7i(z;,A)) =/H,e^,(t;,A), 

(0|J5o(0)|5oW) -/s„, 

(0|J^^(0)|5i(«,A)) =/5,eg,KA). (18) 

As an example, we present the derivation of the sum rule for the coupling constant ff^h^ ^, where p denotes the 
orbital momentum of the final pion. g^hj^ ^ is defined through the decay amplitude for the channel H'l — >• Hi + tt: 

M{H^ ^ ifi + tt) = , (19) 

where q is the momentum of the pion, e* is the polarization vector of the final Hi heavy meson, the isospin factor 
1 = 1, l/-\/2 for the charged and neutral pion, respectively, e^"^ = e'^'^^'^rifj,e*qpVa with e'^'^P'^ the Levi-Civita tensor. 
The correlation function involved in this case is: 

dx e-^'=-(7r(g)|4^(0)j];i(a;)|0) = Re^("^' q.,ViG^^,^Juj,uj') , (20) 

where uj = 2k ■ v and uj' = 2(fc — q) ■ v. When w, lu' <^ 0, we work at the quark level and express the above correlation 
function by the pion light-cone distribution amplitudes: 

G^, (c.,w') = -i f Va e'*(t"+t-') ( (m„ + ma) [2Ai. (a) - Vy (a)] + 2r(a)g ■ 

(21) 

where u = a2+oi^ and u = 1 — u. Furthermore, cj') can be related to g^^h^^^ by the dispersion relation 
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'^H'^H.A^'^) = / ^^^27 ^-w ; -T+ dsi -+ / ds2 - + ---(22) 

Jo Jo [Si - UJ - ie)[S2 - uj' - le) Jq si - uj - it Jq S2 - uj' - le 

with 

P^.^^j5i,s2)-/H./ff.g^,^^^5(si-2Aff.)'5(.s2-2AH) + ... . (23) 
After invoking the double Borel transformation SjJ'^S^?, we extract the double dispersion relation part of Eq. 
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where 
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The function /„ (x) which is introduced while subtracting the contribution of continuum is defined as 

n ^ 

/„(x) = l-e--^^. (27) 



Here we present the details of the continuum subtraction. After invoking the first double Borel transformation to 
the dispersion relation Eq. ([^^ we arrive at 



'^^'^h-hJ^'^')^ ds2e-^^'^^e-^-^^^PH^H,Asi,S2)- 

1 f n f n 



Now the spectral density p{si, S2) can be derived after a second double Borel transformation: 



(28) 



(29) 



According to quark-hadron duality, we can subtract the contribution of the excited states and the continuum and 
arrive at 



-2{toA„h/T-2«oAff/T 



dsi / ds2 e 



-Sicri —S2<72 



BL'^^BZ^B^' B:J G^^.^Juj, uj') , (30) 



where uJc and uj'^ are the continuum thresholds of the mass rules of the H^'' and H doublets, respectively. The terms 

1 1 

of Bu^ B^J G^ijhfj ^(w, w') have general form cu^T'^ = ctr™ / (tri + (72)™^"- Here we assume m, n > to illustrate the 
procedure of the continuum subtraction. 
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(31) 



where s+ = (si+S2)/2, s_ = (s2 — si)/2, 1/T_ = 1/Ti — I/T2 and we assume Wc >^'c- We will work at the symmetry 
point, namely Ti — T2 — 2T and uq = 1/2. This leads to a greatly simplified continuum subtraction: 



T" to! w' 

"77 Tt(2"0 - l)Vri-l+i(7F) 

2™ ^-^ i\(m~ i)\ 1 

<r"/„-i(^), 



(32) 



namely T" r"/„_i(aj^/r). This choice of mq is based on the consideration that the working interval of the Borel 
parameter T of the mass sum rules for H and S is about 0.8 < T < 1.1 GeV[l^, which is very close to that of the the 
mass sum rules for D^{D = H/S/M/T). This will enable us to subtract the continuum contribution cleanly, while 
the asymmetric choice will lead to the very difficult continuum substraction [Tsj . 
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FIG. 4: The sum rule for g^hj^i^ with FIG. 5: The determination of the upper hmit of 

= 2 8 3 3 2 GeV ^ with continuum threshold 

■ ' ■ ' ■ ■ w' = 2.8,3.0,3.2 GeV. 



The binding energy and the overlapping amplitudes extracted in the previous section are used in our numerical 
analysis of the sum rules for the above pionic couplings, together with the following values for A^/s a-nd fn/s UM- 

Ah = 0.50 GeV , fn = 0.25 GeV^/^ , 
As = 1.15 GeV , fs = 0.40 GcV^/^ . 

In this work the tt decay constant is taken to be = 131 MeV. /i^ = m^/(m„ + md) = (1.573 ± 0.174) GeV is given 
in Ref. [l^. The parameters appearing in the tt distribution amplitudes are listed below |1J]. We use the values at 
the scale n = 1 GeV in our calculation. 



0-2 m ^^4 ^00 Woo aio vio hoi hio 

0.25 0.015 -1.5 10 0.2 -3.33 -3.33 5.14 5.25 3.46 7.03 



From the requirement of the stability of the coupling constant to the variation of the Borel parameter T and the 
requirement that the pole contribution is larger than 40%, we get the working interval of T^m < T < T^ax, which is 
plotted in Fig. [S] The resulting sum rule is plotted with = 2.8, 3.0, 3.2 GeV in Fig. S] 

Similarly, we have 
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T/o( 



T 



V2' 



^ r[°l (uo)T/o(^) - 2 {uo) + Vj-'i (^.o) 



■niu + md 
1 



I m . 4~ m J ± II 



m„ + nid 
1 



TO„ + nid 



+ _l!^r["lK)T/o(i^) + Ami Vi-^^o) + V-'\uo) 
■niu + md 1 L II . 
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(33) 



T[GeVJ 



TIGeVJ 



FIG. 6: The sum rule for g^',^ „, with 
continuum threshold cj'. = 2.8, 3.0, 3.2 GeV. 



FIG. 7: The upper limit of T which is not larger 
than its lower limit, indicates the nonexistence of 
a stable working interval of T 
K = 2.8,3.0,3.2 GeV). 



We notice that 
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(34) 



These relations are consistent with the heavy quark flavor-spin symmetry. 

From Fig. [B]and Fig. [7]we can see that there is almost no stable working interval of the Borel parameter T, namely 
the sum rule for g'^hg^ is unstable, so is the case of g'ghjf^- Despite this, we make crude estimates for them, and use 
these estimates to calculate the partial widths of the corresponding decay channels. 



P s p p d d 

i-«17r ^'sfSiTT y^fSiTT 3s'lHi-K yAffHiTT ^'t^Sitt ^'il/f Si tt dT^^HiT! 



p 



0.03 



0.01 



0.7* 



0.7* 



0.2 



0.1 



0.08 



0.1 



TABLE II: The absolute values of the coupling constants. The units of the P- and D-wave coupling constants are 
GeV""'^ and GeV~^, respectively. Here the superscript * indicates the instability of the sum rule for that coupling 

constant. 

The extracted coupling constants and the partial widths obtained with them are collected in Table HIl and Table Hill 
respectively. These numerical values are rather small as a whole. The annihilation of the gluon degree of freedom in 
the decay processes may be responsible for these weak couplings. The tensor structures of the involved decay channels 
are also included in Table IIIIl where the tensor structures of various partial waves are defined as 
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TABLE III: The partial widths of the decay modes D^/B^ -> Z)/B + tt in unit of MeV, together with the tensor 
structures of these decay channels in the heavy quark limit. The masses of the c quark and h quark used in the 

calculation are 1.4 GeV and 4.8 GeV, respectively. 

rnaia2P _ "'^„"10„"2 _|_ "'^„"2/3 ai ^ „aia2 

— 'it "T ^* ' 

rpaia2 „a;i a2 ^ „aia2 -,2 

J-d — It It ~ -^9t It > 

V. CONCLUSION 

We constructed the appropriate interpolating currents for the hybrid mesons containing one heavy quark(q'Qg'). 

Then we calculated the binding energy and the pionic couplings at the leading order of HQET within the framework 
of LCQSR. The mass sum rules and most of the sum rules for the pionic couplings are stable with the variations of 
the Borel parameter and the continuum threshold. For the sum rules for g^hs-jr ^^'^ 9s'^h-k^ ^® ^^"^ "^^^ ^^''^ ^ stable 
working interval of T. We found that the binding energy of the heavy hybrid and M'' are degenerate, so is the 
case of and . As far as the pionic couplings are concerned, the extracted couplings arc rather small as a whole. 

Some possible sources of the errors in oiu- calculation include the inherent inaccuracy of SVZ smn rules and LCQSR: 
the omission of the higher dimensional condensates, in our case the quark gluon mixed condensates etc. in the mass 
sum rules, and the higher twist terms in the OPE near the light-cone, the variation of the binding energy and the 
coupling constant with the continuum threshold ujc and the Borel parameter T in the working interval, the omission 
of the higher conformal partial waves in the light-cone distribution amplitudes of pion, and the uncertainty in the 
parameters that appear in these light-cone distribution amplitudes. The uncertainty in /'s and A's is another source 
of errors to the light-cone sum rules. Finally, the l/mq correction may turn out to be quite large concerning the 
charm quark while such a correction is under control in the case of the bottom quark. 

The weak pionic couplings between the heavy hybrid mesons(gQ5') and the conventional qQ systems render narrow 
partial widths of the corresponding decay channels. We have made a rough estimate of these partial widths. These 
heavy hybrid mesons are found to be quite narrow with a width around several tens MeV in the heavy quark limit. 
We hope that this estimate, together with the calculation on the binding energy of the heavy hybrid doublets, may 
be helpful to the future experimental search of these unconventional heavy mesons. 
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Appendix A: The light-cone distribution amplitudes of the pion 

The 2-particle distribution amplitudes of the tt meson are defined as [T^ 

(0|S(z)7^75d(-z)k-(P)) = iUp^ [\ue'^P'Mn) + y.m^ — Zf, f du e'^^' g^{u) , 

Jo ^ pz Jo 

{0\u{z)^^,d{~z)\7T{P)) = JlII^ f due'^P' 4>p{u) , 

rriu + md Jq 

(0|u(^)a„^75d(-^)k(n) = ipczp-ppzc.) / due'^P'Mu), (Al) 

3 77i„ + rud Jo 

where ^ = 2u — 1, is the leading twist-2 distribution amplitude, (f>[p^cr) are of twist-3. All the above distribution 

amplitudes cj) = {(j^Tr, 4>pi 4>aT9-K} are normalized to unity: jj^ du 4>{u) = 1. 
There is one 3-particle distribution amplitudes of twist-3, defined as [3] 

{0\u{z)a^^j5gsGap{vz)d{-z)\T7^{P)) = i [PaPi^g^p - PaPt^g^p - VfiVp.gva + PjiPuQa^) T{v,pz) , (A2) 
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where we used the following notation for the integral defining the 3-particle distribution amplitude: 

Here a is the set of three momentum fractions aa, and ag. The integration measure is 

j TDa — j daddaudag6{l — q;„ — Ud — ag) ■ 
The 3-particlc distribution amplitudes of twist-4 are 



(A3) 



(A4) 



{0\u{z)jf,j5gsGapivz)d{-z)\'K (P)) = Pf,{paZp - ppZa) — f^rnlA\\iv,pz) + (ppg^^- pagj^,)f^mlAj_iv,pz) , 

pz 

{0\u{z)j^igsGap {vz)d{-z)\TT~ (P)) = pf, {pazp -ppza) — /tt™^ V|| (u, pz) + {ppg^f, ~ Pagpf,)UmlV± {v,pz) , (A5) 

pz 

where Ga/s is the dual field Gap = ^^ajS-ysG''^ . 

We also use the distribution amplitude given in Ref. [l^ : 



Mu) = 6u(l-ii) (1 + 02^2'/' (O) , 

g^{u) = l + (l + ^a2 + 607,3 + y774)C2'/'(0 + (-^02 -67730^3)Cy'(0, 



A{u) 



V±(a) 



(A6) 
(A7) 



^ _ , 16 24 20 
6uu <j — + — 02 + 2O773 + Y »74 



+ 



1 1 

15^16 27 



27 



11 4 \ r^3/2/e\ 



210 



135 



+ ^-y 02 + 21r;4Cj4^ {2u^(10- 15m + Gw^) Inw + 2u^(10 - 15m + Gm^) Inu + mm(2 + 13mm)} , (A8) 



6(m) = ^^(m) - (/!)^(m) , 



0p(u) = I + ( 30r;3 ~ ^ pi ] C',^\0 + { -SV3i^3 - ^ pl ^ ^ P>2 ] C^'^O 



27 



81 



20 ' 



10 

-,3/2, 



-.1/2/ 



(l)^{u) = 6m(1- m) |1 + 1^5773 - ^773^3 ~ - ^P^a2 ) C2'"(0 

T(a) = SeOrisa^adal j 1 + W3 i (7ag - 3) I , 



V||(a) = 120a„adag(uoo + vioi^ag - 1)), 
^^ii (a) = 120a„Q!(iagaio(Qd — a„) , 



-30a^ 



-4x(q;) = 30a (a„ - ttd) 



/ioo(l - ag) + hoi ag{l - ag) - Ga^ad 
1 



"9(1 -"s)- +arf) 



/ioo + ^-01% + 2 ^io(5ag - 3) 



(A9) 
(AlO) 

(All) 

(A12) 

(A13) 
(A14) 

(A15) 



where C™(^) are Gegenbauer polynomials. 

The definitions and the specific forms of the rj light-cone distribution amplitudes adopted in the text are similar to 
those of the pion. For more details see Ref. [l^ . 



